Gravitational self-force in the ultra-relativistic limit: The "large-N" expansion 
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We study the gravitational self-force using the effective field theory formalism. We show that in 
the ultra-relativistic limit 7 — > 00, with 7 the boost factor, many simplifications arise. Drawing 
parallels with the large N limit in quantum field theory, we introduce the parameter 1/N = I/7 2 
and show that the effective action admits a well defined expansion in powers of A = Nt at each 
order in 1/N, where e = E m /M with E m — 7m the (kinetic) energy of the small mass. Moreover, 
we show that diagrams with nonlinear bulk interactions first enter at 0(1/ N 2 ) and only diagrams 
with nonlinearities in the worldline couplings, which are significantly easier to compute, survive in 
the large JV/ultra-relativistic limit. As an example we derive the self-force to 0(X 4 /N) and provide 
expressions for conservative quantities for circular orbits. 



Introduction. During the last years a new formalism has 
emerged, based on effective field theory (EFT) ideas bor- 
rowed from particle physics, to study binary systems in 
General Relativity. Originally proposed to tackle the 
Post-Newtonian approximation for non-spinning [1] and 
spinning [2] inspirals, the EFT approach has produced 
a breakthrough in our understanding of gravitationally 
interacting extended objects within a relatively short pe- 
riod of time [TMl9|. Meanwhile, EFT ideas were also ap- 
plied (besides particle physics) to different areas, such as 
cosmology [20H22] . electrodynamics [23], fluid dynamics 
[2"4l |2"5] , and in particular to extreme mass ratio inspirals 
(EMRIs) [2l)H2"g] . which is the subject of this letter. 

The study of the self-force problem within the EFT 
approach was initiated in [26] where power-counting and 
leading order effects were worked out and a proof of the 
effacement of internal structure for EMRIs was given. As 
it is traditional in the EMRI community, perturbative 
calculations are performed in powers of q = m/M where 
m represents a small mass object orbiting a much larger 
black hole with mass M. To date, only second-order 
0(q 2 ) equations of motion are known [29| [30] . 

In this letter we study the ultra-relativistic limit, i.e. 
7 — > 00, of the self- force problem. Inspired by an analogy 
with the large N limit in quantum field theory |31j . we 
show that many simplifications arise that are not cap- 
tured by the (traditional) m/M power-counting. Wc 
show that, introducing the expansion parameter 1/N = 
7~ 2 and defining A = Ne with e = E m /M and E m = 7m, 
the gravitational effective action (which yields the self- 
force) admits an expansion of the type 

See = L/N(l + \ + A 2 + ...)+ 0(\/N 2 ), (1) 

where L ~ E m M(— jmM) is the angular momentum of 
the small mass (in Gn = c = 1 units used throughout). 
A similar expansion applies to the one-point function, 
h fJil/ (x), which can also be used to compute the self- force. 
Our goals here are: to derive the new power-counting 
rules in the large N limit, to show that diagrams with 
nonlinear bulk interactions are subleading in the 1/iV ex- 



pansion, to report the gravitational self-force to fourth 
order in A at leading order in 1/N, and to provide expres- 
sions for conservative quantities for the case of circular 
orbits. We conclude on a more formal note with some 
comments on the problem of finding the self-force in the 
exact massless limit for a photon moving in a black hole 
spacetime. 

Power- counting rules. Our setup is the same as in 
the standard EMRI case except that we consider ultra- 
relativistic motion where the boost factor 7 is large, 

7 = l/yj-g^v" > 1. (2) 

Here, v a = dz a /dt, z a is the small mass' worldline co- 
ordinates, t is the coordinate time, and g^ v is the back- 
ground metric of the black hole with mass M. In order to 
achieve this condition we imagine the mass m in a bound 
orbit near the light ring in Schwarzschild spacetime or 
moving in a (very) rapid fly-by. 

We next find the scaling rules of various leading or- 
der quantities. The orbital frequency is related to the 
wavelength of the gravitational radiation through CL> or b = 
dcj)/dt ~ 1/A gw . Because A gw ~ A I it follows that dt ~ M 
(also dx % ~ M). Hence, the proper time along the ob- 
ject's worldline scales like dr ~ dt/-f ~ M/7 and its 
four- velocity is u a = dz a /dr — jv a ~ 7, for an ultra- 
relativistic motion. For the scaling of the metric pertur- 
bations h^v produced by this ultra-relativistic small mass 
m we use the leading order solution 

VW~ / G ava , p ,(x,x')T^(x'), (3) 

J x' 

with T a/3 (x) cx mj drS 4 ^ - z' j (t))m q u' 3 / % /=^, and 
/ = J d 4 x^/^g. We find hp V ~ E m /M = e where we 
used V Q ~ d a ~ 1/M and G^a'p* ~ 1/M 2 for the scal- 
ing of the Green function in a curved background (this 
follows almost entirely from dimensional analysis). Fi- 
nally, the leading order effective action scales like 

0*1 =-m I 'dr ~ mM/j - L/N, (4) 
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as anticipated. The scaling rules are summarized below: 
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Because of these rules, the condition that perturbation 
theory is under control in the ultra-relativistic limit de- 
mands not only e to be small but also 7 2 e <C 1. The 
reason is simple. After including the perturbation the 
point particle action scales like 



S 



PPL 



(It 



^Jl-^h^v^, (5) 



where we used According to the scaling rules, 

hfj,i> ~ e : hence we must require A = ej 2 = eN <C 1 
for the perturbation to be considered small with respect 
to the background. This is the regime of validity of our 
approximations. In other words, we formally take the 
limit e — > 0, N — > oo with A fixed and small. 

To obtain the different scalings for the possible terms 
that contribute to the self-force we first need to isolate 
the building blocks of our Feynman diagrams and power- 
count them. We have either worldline or bulk vertices, 
which we summarize next. Using our power counting 
rules we have for the vertex describing the interaction of 
the small object m with n gravitational perturbations: 

*"v* - M (^) (7 2 )" = mM^-\ (6) 
which arises from expanding the point particle action 

oo „ 

S PP = m ^2 c n / dr (/i Q(9 u q m /3 )™ , (7) 
n=o J 

where the {c„} are dimensionless numbers. Notice we 
truncate the external legs and we do not yet include the 
scaling for h a p , which ought to be contracted with world- 
line or bulk couplings and will introduce an extra factor of 
G^ua'f}' ~ M~ 2 for each propagator in a given diagram. 
Next, we need the bulk vertices that follow from expand- 
ing the (gauge-fixed) Einstein-Hilbert action in powers of 
hfj, v about the given background spacetime g^ v . At the 



,th 



order this is given schematically by 



OO r. 

5 eh = / VhVhh n ~ 2 , 

n=2 J x 



(8) 



where V indicates the covariant derivative. It is easy 
to show that the vertex for n interacting gravitational 
perturbations scales as: 

2 n - 1 
1 n 

in four spacetime dimensions. This completes the power- 



counting rules for the building blocks of the EFT for- 
malism. To compute the classical effective action we 
simply need to add up all possible tree-level diagrams. 
(By this we mean we do not include closed gravitational 
loops that represent quantum effects.) The effective ac- 
tion then takes the form 



S eS [z" 



Using the rules previously derived we can power-count 
each diagram in the effective action, hence their contribu- 
tion to the self-force. We show next that only diagrams 
without bulk nonlinear interactions survive in the large 
TV limit. For that purpose it is illustrative to compare the 
scaling of the following diagrams, which enter to C(A 3 ): 
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(10) 

(11) 
(12) 

(13) 



We already start to see the pattern: bulk nonlinearites 
are suppressed in the large N limit. For a generic con- 
tribution let us consider a diagram with N m mass inser- 
tions, bulk vertices with fc-legs and N p propagators 
(including internal ones). From our power counting rules 
we obtain the scaling 



(Mm//» M 2 « 



-N p ) 2(2N p -J2 k kV k ) 

I 1 



(14) 



where iV* ot = ^ fc is the total number of bulk vertices. 
Let us first look at diagrams with 2V* = 0. Using: N m + 
iV* ot — N p — 1 = 0, which follows from the topology of 
the diagrams that contribute in the classical limit, the 
expression in ( 14 ) turns into 



L 



-1) 7 2(JV„ 



D = ii A (^i) (^ = 0), (15) 



N 



and is thus a 1/N contribution. From a given order in 
A (namely, N m fixed) adding bulk vertices (and internal 
propagators) will only introduce powers of 1/N (see ( 14)) 
since we need at least two bulk vertices to increase the 
number of internal propagators and each bulk vertex has 
at least three legs. Intuitively this is simply because, for a 
fixed number of mass insertions, we lose powers of N from 
propagators attached to two worldline couplings, which 
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are promoted to a bulk interaction. This is transparent to isolate the part of (17) that produces the divergences 



in the terms depicted in ( 10 )-( 13 1 



The gravitational self-force in the large N limit. Self- 
force effects in EMRIs are intrinsically non-local, de- 
pending on the past history of the small object's motion 
around the larger black hole. Capturing these real-time 
dissipative interactions with an (effective) action requires 
a careful handling of Hamilton's variational principle of 
stationary action so that it is consistent with initial value 
data for open system dynamics (i.e., the motion of the 
small mass). This issue was emphasized in [32] where it 
was motivated by the classical limit of the "in-in" for- 
malism (33) . A rigorous framework to handle this in a 
completely classical context was developed in [34] and 
applied to derive radiation reaction forces through 3.5 
post-Newtonian order using the EFT method in [S]. We 
will elaborate on the details of this construction for the 
self-force problem in the large N limit elsewhere [3"5] , 

As we have shown, in the ultra-relativistic limit we 
can ignore all self-interactions of the metric perturbation 
that do not happen on the worldline. This means that 
the action for the small mass object and the metric per- 
turbations can be taken as 



S[z",V] 



1 

64?r 



At \ 1 — h a pu a uP, 



(16) 



where we imposed the Lorenz gauge for trace-reversed 
perturbations. For the reader worrying about finite size 
effects (e.g., Ce f At £ a p£ a ^) one can easily show they 
are highly suppressed in the large N expansion, first en- 
tering at 0(X 4 L/N 5 ). This has important consequences 
in the regularization of the theory because, as we shall 
argue, we will not encounter logarithmic divergences but 
only power-law, which will be handled via dimensional 
regularization (and set to zero since they involve scale- 
less integrals). We briefly discuss below the general pro- 
cedure for calculating the relevant diagrams in the ultra- 
relativistic limit, and we will elaborate on technical issues 
elsewhere [35 . However, these steps follow closely those 
taken in a nonlinear scalar field model of EMRIs 28 . 

Computing the surviving diagrams in the effective ac- 
tion, or the diagrams for the one-point function h^ v {x) 
below, involves worldline integrals over the retarded 
propagator, 

/CO 
At"G%, y , s „(z»\z»")u"<"u s ", (17) 
-00 

which are in general divergent. Here, a prime on an index 
indicates the point or proper time that the quantity is 
being evaluated so that u a — w Q (r'), u 1 = u 7 (t"), 
etc. Following 36j we split this expression into a regular 
^V/3'7"<5" an d singular G^,g, „ s „ piece, which allows us 



It is useful to write the singular integrals in a momentum 
space representation, which can be given whenever the 
two points on the worldline can be connected by a unique 



geodesic. Using the above decomposition one writes (17) 



as Xs(z M ) +/_r(z m ) where the singular and regular parts 
are, respectively, given by 



Is(z^) = 4w Q ^'P Q ^ vs -(^')Re 

f°° d d k e -^V"-T') 

x 



At" ujl'ufi' 



_ x (2ir) d (fc0)2 _ p + 



(18) 



ic 



and 



I R {zP ) = u a u^ I At" D$^, s „ (*" , )uV (19) 



6" 



where 



= 6(r" - r out )9(r in - /')G^^,«(^', *"") 
+ e(r out - t")Q(t" - 7i n )G^ v ,,„(^\ *""). (20) 



The integral in ( 18 ) is written in A spacetime dimen- 



sions in momentum space where the momenta are dual to 
Fermi normal coordinates, and u7, = g 1 \" (^ > z ^ 
is the result of parallel propagating the velocity vector at 
z m (t") to z^(t') using the propagator of parallel trans- 
port g 1 \"{z^ ,z M ). Also, Ti n (and r out ) are the proper 
time values at which the worldline enters (and leaves) 
the normal neighborhood of z^(t'). See [37] for more 
details about bi-tensor calculus. As we mentioned, the 



singular term in ( 18 ) is easily shown to vanish in dimen- 



sional regularization because it is a (scale-independent) 
power-law divergent integral. As a consequence, the reg- 
ularization of the theory becomes straightforward in the 
large N limit. 

As outlined in 28J , in a theory that has only worldline 
interactions as the relevant couplings, it is simpler 
to compute the metric perturbations at a field point, 
hfj,v{x), rather than the effective action since we can 
simply substitute the resulting perturbative expression 
into the worldline equations of motion (renormalizing 
parameters if necessary) to compute the self-force. We 
show the results next. 

Metric perturbation & self-force to 0(X i /N). In the 
ultra-relativistic limit the diagrams contributing to the 
one-point function can be represented as the convolution 
with a worldline master source, 

M*) = / AT'G% a , ,(x,z»')S^'(z»'), (21) 

that is completely finite and given through 0(X 4 L/N) by 
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computing the diagrams 



from which we find the master source to be: 
5«V) = y^'{l + f/fl(^') 



l 2 

3m 3 



128 7 7 6 " J " 7 128 

3m 3 i i if ti a a 

+ -^IR^ )^U S jdT"D R , s , e „ r ^vP I R (z» ) 

m 
+ 64 



u<u 5 'l dT''DX, s , t ,, v ,,u e ''u r, " u p" u X' 



dr"'D R „ a „,vJ v° I R {z» ) + 



(22) 



From the master source we can then compute the regular 
part of the metric perturbation evaluated on the world- 



line, i.e. h*(z^), by convolving (22) with D 



(20 1. In order to obtain the self- force we compute the 



worldline equations of motion via the variational prin- 



ciple on (16 1. We thus arrive at the desired result, in 



principle valid to all orders 
g^a + h R , 



i R 



1 - h R s u~tu s 



2 1-ft £ >V U )Vuh aft , (23) 

where h R v is evaluated on the worldline using the master 
source in ( 22 ) , and we have absorbed a divergent piece 



into the mass m. (These divergences are set to zero in 
dimensional regularization; recall that there are no other 
counter terms at leading order in l/N.) The formal per- 
turbative expression for the self-force can be easily found 



by expanding out ( 23 1 to the desired order (modulo some 
ambiguities discussed in detail in |35j). 

Conservative self-force for circular orbits near light ring. 
On circular orbits in a Schwarzschild background the 
symmetry of the system implies that the regular integral 



in (19) is a function only of the orbital radius r Q , and is 



independent of time as long as we consider only the con- 
servative part of the self-force, implying I^z^) = Ir(t ). 
The master source in (22) then becomes, to 0(X 4 /N), 



' N 



rn 

~2 



Ir 



5m 3 
~6T 



1 r. 
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(24) 



We can now use the above master source, and corre- 
sponding regularized metric perturbation h R v evaluated 
on the worldline, to obtain an (implicit) equation for the 
energy of the orbit, E, by solving 



1 



3M 



2M 



,E Z = 



(fr,V) 2 

7-2(1 - hf € u s u e ) 



(25) 



where we used the radial component of the (non- 
perturbative) equation of motion in ( 23 ) . (The u a and 



h R a depend implicitly on E.) In principle we need to ex- 



pand ( 25 ) perturbatively in powers of A about the back- 
ground energy Eq of a circular geodesic. This will be 
given in detail in |35j . Unfortunately, we find that there 
are contributions (from h R r (r ) and h R a (r a )u a ), start- 
ing already at 0(X 2 /N), which have not yet been ob- 
tained numerically in the literature. To this extent, we 
expect our results will encourage the community to com- 
pute these terms in the future. 

Concluding remarks. We have introduced the large N ex- 
pansion for computing the gravitational self-force in the 
ultra-relativistic limit and shown that, at leading order 
in l/N, it reduces to a (mostly) combinatorial problem. 
As an example, we derived the self-force through fourth 
order and gave the (non-perturbative, implicit) expres- 
sions for the conserved energy for circular orbits near the 
Schwarzschild light ring. Our results are most useful the 
larger 7 is, provided Ne = j 3 q remains fixed and small. 
For example, in our computations, ignoring l/N 2 correc- 
tions requires X 2 /N 2 < A 4 /TV, or I/7 4 < q, while at the 
same time j 3 q < 1 for perturbation theory to stay under 
control. Therefore, our range of validity lies somewhere 
between I/7 4 < m/M < I/7 3 . This window obviously 
increases the less accuracy we demand. 

The gravitational self-force has received significant at- 
tention lately, due to some surprising agreements with 
numerical results outside its range of validity (formally 
replacing m/M -> mM / (m + M) 2 ) [38-40 . These com- 
parisons, however, only relied on leading order self-force 
effects. Our results in this letter open the door to check 
and improve such computations to very high orders in 
the large N limit. As it is often the case, these approx- 
imations may shed light on the dynamics in scenarios 
where 7 is not significantly large, and perhaps even in 
cases where the mass ratio is not taken to be small. We 
leave this road open for future work. Our results should 
also be useful to further calibrate semi-analytic merger 
models from the ultra-relativistic regime (e.g., see |4"T]). 
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Let us finish by commenting on a more formal aspect of 
the ultra-relativistic limit. As it is well known, a boosted 
Schwarzschild black hole turns into an Aichelburg-Sexl 
(AS) Shockwave in the ultra-relativistic limit with e fi- 
nite |42j . One simple way to recover this solution is 



using Polyakov's action (Sp i y = J da{x(a) 2 /e 



! }) 



[43] . which is finite in the massless limit. A special fea- 
ture of this point particle action is that it does not in- 
troduce worldline non-linearities, only bulk-type which 
are present through the Einstein-Hilbert action. How- 
ever, all the non-linear terms cancel out for the AS so- 
lution, which is linear in Gm [HI- This is not the case 
in a black hole background (with finite mass M) because 
the Shockwave can encounter its own "echoes" [44]. It 
would be interesting to study AS Shockwave dynamics 
in non-trivial backgrounds as another approach to the 
ultra-relativistic self-force, for instance, to study the dy- 
namics of light crossing a black hole, the merger process 
in binary systems |45j or to understand high-energy grav- 
itational collisions [3B1I17]- (For the case of photons, it 
would also be instructive to compare with the geometric- 
optics limit of the Einstein-Maxwell equations.) While 
this is not the same limit studied here, it would be in- 
teresting to understand the seemingly dual relationship 
between both approaches and the connections (if any) 
between worldline and bulk non-linearities. 
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